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Abstract 

We show that the Bel-Robinson (BR) tensor is - generically, as well as in its original GR 
setting - an autonomously conserved part of the, manifestly conserved, double gradient of a 
system’s stress-tensor. This suggests its natural extension from GR to matter models, first to 
(known) massless scalars and vectors, then to massive ones, including tensors. These massive 
versions are to be expected, given that they arise upon KK reduction of massless D -\- 1 ones. 
We exhibit the resulting spin (0,1,2) “massive” BR. 


1 Introduction 


The motivation for the present Bel-Robinson (BR) [1] exercise is twofold: First, to dispel any 
lingering mystery about these conserved four-index tensors’ existence by showing that one of their 
avatars is just a - separately conserved - part of a stress-tensor’s double gradient, exactly as is 
the case in GR [2]. Second, to exploit this fact by extending BR first to massless, then to massive 
matter. Our explicit examples will include scalars, vectors (Proca) and massive tensor - Fierz- 
Pauli (FP) - spin 2, as well as a “how-to” outline for nonlinear - massive gravity’s (mGR) BR. 
The existence of massive BR extensions is separately suggested by the fact that massless systems 
in D -|- 1 become massive in D under a periodic KK compactification of the “fifth” dimension [3], 
even though, as we will see, the D = 5 BR neither does nor should ~ unlike the action or the field 
equations - reduce directly to the massive D = A form. 
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2 BR-stress tensor connection 


Consider for concreteness the free, initially m = 
conserved is 

(p^y 


0, scalar field in any D, whose on-shell (□(/> = 0) 


2 *^,7 ‘/’7 ’ 


( 1 ) 


obviously da 9/3 T^y remains conserved on {nu). Slightly less obvious, but fundamental to our aim, 
is the separate conservation of its “spread-out” derivatives and the remainder, defined by 


9q 9/3 T^v — 4^,ua “f (p,^a 4^,u/3 'P'^a 4^,1/ “f (p,^i 4’,vf}a 4’,'yfici > 

= + [^^lua|3\ 

Their independent conservation is trivially verified explicitly, or by noting that A’s higher derivative 
structure requires it. We will see that the same holds for higher spins, so henceforth we drop A, 
and concentrate (mostly) on the “spread-out” parts of 99T. 


The extension to massive scalars follows the same prescription, namely adding the “spread- 
out” da 9/3 acting on the mass part of here —Xj^vr? rj^ycp'^ in our mostly plus signature, where 
(d—m?) (p = 0. The compactification of with periodic m-dependence of the “5th” dimension 

is of course precisely (Dd — rrp). For completeness, we display the massive 

B^yaP — 4*,va T (p,lice *^,7/3 4’’^a ^ 4^,a 4^,13 (3) 

We may now explain the disparity mentioned above between the D-components of and 

the massive i?((7o/3 ™ That they should not agree is obvious from the fact that their conservation 
properties differ, since only 

= 0, (4) 

and of course the second term does not vanish. It is easy to verify that this “discrepancy” between 
these two BR is precisely what is required for their respective conservations. 


3 BR for s = 1, Proca fields 


The higher spin equivalents of our scalar system follow exactly the latter’s pattern above. For 
vectors, start with = F^aFy'^ — X/Arj^j^y Fap F°‘^, apply the 9 q, 9/3 prescription, and add the 
da 9/3 of AT^jy = w? (Af^ Ay — X/2rifj_y Aa R"). The conserved (on {pLv) only) form of the Proca BR 
is thus 


F^yapi-^p) — F^a,a Fy 


FF^a^p Fj^ Q, ^ rj^y 


: FfjY^pFfTi (^A^ a T ^u,a -^cr,l3 -F a) ■ 


(5) 

For m = 0, it agrees with the original Maxwell version [4]. Conservation of B in all cases (including 
GR) consists of two separate parts: one uses the “Bianchi” identities, here the cyclic 9[^ Fy^^ = 0, 
while the other depends on the dynamics, here 9^ F^y — m? Ay = 0, and similarly for the scalars 
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treated above. The mass terms only affect the dynamical part, of course. 


4 Spin 2 


Our new example is the massive spin 2 - Fierz-Pauli (FP) field, consisting of linearized GR as 
kinetic, plus a - unique - mass, term. This model has received enormous recent attention as the 
starting-point for so-called massive gravity (mGR), its nonlinear extension with full GR plus a 
non-derivative mass term, necessarily involving a fixed background metric. Rather than proceeding 
via the (linearized) BR route, it is instructive to first start with the FP stress-tensor, as for lower 
spins. The latter, unlike that of linearized GR alone, is perfectly well-defined, there no longer being 
any gauge invariance to respect. We will use a - legitimate - shortcut below, namely since we are 
only interested in and its related ddT “BR” forms on-shell - they are not dynamical currents 
- we apply the five constraints = 0 = on the h^y field that reduce it to describe just the 

2s -|- 1 = 5 spin 2 excitations. As a result, the action is simply 

1(h) = - 1/2 j d‘^x[(hf,y^a)‘^ + vf? (h^y)% ( 6 ) 

here all indices are moved by the - mostly plus - Minkowski metric so the field equations 
are (□ — w?)h^y = 0. To obtain the stress tensor, we could follow the Belinfante procedure; turn 
all rj^y —)■ g^y(x)^ vary the action with respect to g^^y, then set them back to to define T^j^y. 
In addition to the g^'^ and the overall ^/—g factor in ( 6 ), a new (higher spin) phenomenon arises 
because the partial derivatives da now become covariant ones. Da = (da — F) acting on the tensors 
h^y. Their variations, AT ~ d(Ag) give rise to novel contributions to T^^y, namely ~ d(hdh), that 
is not only ~ dhdh, but also ~ hddh, absent in the lower spins. A dull calculation yields 


Tfiu — T(g)^y -\- T{T)^y 

'dd(g)^y ‘^h^a,l3hya,l3 T hafj^^ ha/3,y X hlii'(^ci h/S'f) T 'IXl 


2 h^a hya „ ( ^"fS ) 


(7) 


T(Y^)^y hj 3 ^(h^^,yi 3 -\- V -f-)- g) (hy^,i 3 hg^,^ -t- zz -f-)- fx) ^hy^^gh^^y,g (hy^\Z\h^'^ -t- g ix) 


which is easily verified to be conserved on-shell, namely with (□ — m?)h^y = 0 = d^ h^’^ = 0 = 
/i“. This T^y is of course not unique: there is always an ambiguity - for all stress-tensors ~ 
of adding identically conserved “superpotentials” dj 3 dsH'^^^^'^’''^\ that also do not contribute to 
the integrated (Pfj.,Jfj.y) (and of asymmetric of the type present in the Einstein pseudo¬ 

tensor). Since we are concerned with the “spread” dad/s of T^^y, the hddh terms present a potential 
obstacle. Instead of removing them by finding a suitable superpotential, we use the already available 
Landau-Lifshits (LL) pseudotensor, since it embodies all the needed properties, being symmetric, 
conserved on-shell (also for m 7 ^ 0 , as can easily be checked from the detailed form below) and 
- most important - of the desired dhdh form. Including the mass contribution, and using hf^y's 
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traceless-divergencelessness, it reduces to: 

2 hpp cr hi/j^ p hjjpp hpjj p hj^p i/ + hpp^p hyp p -\- — hp\ p hpx ^, — 'f]pv{hp\,i3) 

+ m \hpa hi/a ~ ^ Vpi' ^0/3] ~ ■ 

( 8 ) 

Now we may safely proceed to take its “spread”, r^(LL), 4-index extension as for 

lower spins. We need not spell out the obvious result, ~ 9Q9^r^(LL) ~ (ddh) (ddh) + 

m? {dhdh), with the assurance that it is indeed conserved on FP shell (as is easily checked as well). 


The alternate route to the above “low spin” procedure is to return to FP’s GR roots and use 
a true (linearized, of course) BR kinetic term, namely: 
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Bpua/3 — Rcrfira BuutP T Bapirfi Bcrura - g gp-v OaP RxarS, 
whose divergence is 


B 


0 


^(TVTp (-^crT,a “1“ ^ ^ P 


using the Bianchi identity. [This form of BR is equivalent to 

^ -^cri^r/3 “1“ R(t^t(3 ^(jutol RcxarS ^/SarS 5 

(only) in R = 4 where Rpa/S'^ RuaP'y = j 9pu Ra/S^S Ra/S^s]- On FP shell, 
Rpua0,p — 2 ^ {hl3i/,a ^ 01 /,/s) Ruf) — 2 ^ 

we get 

Rpual3,p 2 ^ Rcri/T/S (hcFT,a ^(70,t) T Ck /?• 

Now we add to B^ of (9) or (11) the mass term 




gpv {h CXCF^T hpa,T - h crr,Q hpa^r - h CTT^^ ^Q;(T,t) 


(^Q;/i,(T ^/3z/,cr ^Q;/i,(T ^cri^,/? ^/3z/,cr T hp(j,a “1“ OL ■f-)' /5) 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


(14) 


It is then easy to prove that B^=B^ + AB is conserved on {piv). 

The second massive model is, as we mentioned, massive gravity (mGR); it adds to the full 
GR action a 3-parameter family of mass terms that depend on both the full “metric” g^y of the 
Einstein part and on a fixed background fpy. This hybrid is NOT of course a Riemann space theory 
but a (complicated) spin-2 field action in a fixed, say flat, space fpy, where local stress tensors are 
well-defined. What changes from (9,11), then, is that is now the full nonlinear one, as is the 
AT ~ contribution. The dependence on being extremely complicated for all three allowed 
(by ghost-freedom) mass terms, we have not bothered to compute their AT^y - they are sufficiently 
sick [5] to warrant omission. Nor have we considered the last variant in which both (g^y and fpy) 
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are dynamical, their action being the sum of their full Einstein contribntions pins a non-derivative 
conpling term ~ Again, since this model is likely to share the acausality of mGR, we 

have not looked for its BR. These formal tasks we are happy to leave open. 


5 Summary 


We have demystified and generalized the venerable BR tensor beyond its GR roots, both by express¬ 
ing its (known) avatars for massless scalar and vector matter fields as conserved double gradient 
parts of their stress tensors, and by generalizing it further to include their massive extensions for 
spins (0,1), as supported by KK considerations. For spin 2, we exhibited free massive spin 2 FP’s 
BR and noted how it could be obtained for nonlinear massive gravity models. Onr process was 
greatly simplified throughout by the separate conservation of the two disparate parts of ddTfj,^, 
allowing us to concentrate on the more physical, “spread derivatives” component. 

Finally, while we have worked in flat space for simplicity, snitable covariantizing of the mass¬ 
less scalar and vector BR is known.The resnlting total BR ~ with extra non-minimal terms required 
by non-commutativity of covariant derivatives - are now (covariantly) conserved on combined Ein- 
stein-l-matter shell. The mass terms’ contributions are now of conrse to be added as ~ D D AT^^. 
We have not carried out FP’s (or even linearized spin 2) BR covariantization, bnt it may present 
a new set of problems, since at least massless tensor fields cannot consistently couple [6] to GR, 
except of course when is the linearized metric deviation about the background! 
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